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A STUDY OF THE FEASIBILITY OF APPLYING CAPACITIVE 


DISPLACEMENT-MEASURING TECHNIQUES TO 
OPEN-MESH GRID STRUCTURES 

By Richard DeLoach 
Langley Research Center 

SUMMARY 

This paper provides theoretical and experimental support for the hypothesis that a 
noncontacting capacitive displacement-measuring transducer can be used to measure dis- 
placements in open-mesh grid structures with extremely small surface areas. Examples 
of such structures include reflector nets which form the primary substructures of radio 
telescopes and radar receiving units. Test models of this kind of structure are typically 
so lightweight and flexible that mounted transducers would alter their dynamic response 
characteristics. Noncontacting transducers are therefore required for modal analysis 
and structural stability tests on this class of structure. 

The purpose of the paper is to evaluate the feasibility of employing a capacitive 
probe in displacement measurements of structures with low surface areas, particularly 
open-mesh grid structures. A solution is found for the field of a charged square-mesh 
grid mounted parallel to a flat grounded plate. A variational technique is used to deter- 
mine the capacitance of this grid-plate capacitor in terms of easily measured geometric 
parameters. It is known that the capacitance of a parallel-plate capacitor is directly 
proportional to the area of the plates. Because the relative surface area of a large-mesh 
grid is extremely small, original estimates of the capacitance of a grid-plate system were 
correspondingly small. However, the theoretical capacitance formula derived herein 
predicts grid-plate capacitance values which are orders of magnitude larger than these 
preliminary estimates. The theory is supported with experimental data. Results show 
that a noncontacting capacitive transducer is indeed feasible for detecting displacements 
in surfaces with extremely small surface areas. 

INTRODUCTION 

The purpose of this report is to demonstrate that a noncdntacting capacitive 
displacement-measuring transducer can be used to measure displacement amplitudes in 
open-mesh grid structures possessing minute surface areas. Examples of such 
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structures include the reflector nets which form the primary substructures of radio 
telescopes and radar receiving units. 

In July 1969 a need arose at the Langley Research Center for an instrument to mea- 
sure the dynamic response of two models of a radio telescope satellite called LOFT (Low 
Frequency Telescope). Both models have essentially the same overall configuration, 
which consists of a spin- stabilized reflector net fabricated from an open-mesh grid of 
flexible, very lightweight strands of rubber-coated steel yarn. The larger model has a 
diameter of 15 meters while the smaller model has a diameter of 5 meters. The pro- 
posed testing program calls for the models to be rapidly rotated in a large vacuum cham- 
ber with the resulting centrifugal force used to deploy the reflector net and maintain its 
shape. (It is necessary to conduct the tests under vacuum to avoid the problems of air 
drag in deploying the model.) Bursts of forced air from jets mounted on independently 
rotating exciter arms under the reflector net will be used to excite resonance modes in 
the net. It is these mode shapes which are to be measured. The tests are complicated 
by a number of factors. Among the most serious complications are the very small sur- 
face area of the reflector net and the fact that it is extremely flexible and lightweight. 
Conventional mounted transducers are not practical since even the smallest ones avail- 
able are massive enough to alter the dynamic response of the structure. 

It has been suggested that a small, flat plate might be mounted near the reflector so 
that the reflector net and plate will form a capacitor. As the grid-plate geometry changes 
as a result of the modal displacements of the net, the grid-plate capacitance will also 
vary, and these changes in capacitance can be related to the motions of the grid. 

Capacitive probes have been in use since 1920 when R. Whiddington and his collab- 
orators first investigated their potential application to the measurement of pressure, 
temperature, and force, as well as displacement (refs. 1 to 4). Since then, capacitive 
probes have found an ever-increasing range of applications (ref. 5). However, most of 
the standard capacitive measurement systems in use today involve very close proximity 
between the probe and the subject surface (refs. 6 and 7). This close proximity causes a 
maximum percentage change in capacitance for a given small displacement in the subject 
surface, providing a high degree of resolution in the measurement system. In the LOFT 
application, however, close proximity between probe and subject is not practical. Dis- 
placement amplitudes on the order of 5 to 10 cm are expected in the tests, which means 
that the probe must be mounted at least that far from the subject. Actually, in order to 
avoid fouling the net during model deployment, the probes should be at least 25 cm from 
the surface of the net. Most capacitive displacement-measuring systems involve probe- 
subject separations on the order of millimeters or less. Not only are the probe-subject 
separations necessarily large in the LOFT application, but the relative surface area of 
the steel yarn which makes up the net is minute. (Each strand is only 0.25 mm in 
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diameter, while adjacent strands can be on the order of 150 mm apart.) For a given 
probe, the probe-subject capacitance is proportional to the area of the subject and 
inversely proportional to the probe-subject separation. Because of the large probe- 
subject separations and minute subject surface area involved in the test, it was originally 
feared that the available probe-subject capacitances would be much too small to be prac- 
tical. It was therefore necessary to determine analytically how the grid-plate capaci- 
tance would depend on grid-plate geometry parameters in order to evaluate the feasibility 
of employing a capacitive transducer in the LOFT tests. 

A square-mesh parallel grid-plate capacitor is used as a model for the analysis. 
The field of such a capacitor is solved, and a variational technique is used to determine 
the capacitance of the grid-plate capacitor model from the field solution. The theoretical 
results are compared with experimental data. 

SYMBOLS 

A n capacitance geometry constant, Dn/V 

a mesh parameter, distance between adjacent parallel wires in grid 

a n general Fourier coefficient 

b grid-plate separation 

C capacitance 

C t capacitance per unit length of a parallel wire-plate capacitor 

D n potential distribution constant 

d n radial distance from an arbitrary point in the plane of a grid of parallel 

charged wires to the axis of charge of the nth wire 

d n ' radial distance from an arbitrary point in the plane of a grid of parallel 

charged wires to the axis of charge of the nth image wire 

L length of a rectangular parallel-surface capacitor 

n positive integer 
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q excess charge 

r wire radius 

rp A 

_ 8neV 

U electric field energy 

V potential difference 

W width of a rectangular parallel-surface capacitor 

x,y,z Cartesian space coordinates 

a separation-of-variables constant 

6 displacement of the axis of charge of a charged wire due to the influence of 

an adjacent grounded plane (see fig. 4) 

e electric permittivity 

A linear charge density 

r volume of integration 

c p potential distribution 


ANALYSIS 

Figure 1 is a sketch of the grid-plate capacitor which is used as a model for the 
analysis. The grid is a square mesh of cylindrical wires, each of radius r, separated 
axis-to-axis by a distance a. The grid is parallel to the plate and is separated from it 
by a distance b. The grid and the plate have the same length L and the same width W. 
In the analysis, the field of this grid-plate capacitor is solved and a variational technique 
is used to determine the capacitance from the field solution. This variational technique 
is described in the following section. 
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Approach 

To invoke basic definitions in determining the capacitance of a pair of arbitrary 
conductors held at some fixed potential difference, it must be possible to determine the 
total excess charge residing on either of the conductors. For systems possessing rela- 
tively complex geometries, this can be a rather difficult task. However, there do exist a 
number of techniques for determining the potential distribution associated with a charged 
system. It would therefore be useful to be able to express the capacitance of some arbi- 
trary charged system in terms of the potential distribution of that system. This can be 
done by equating the energy stored in the field of an arbitrary capacitor to the energy of 
an arbitrary electric field. If U is the field energy, C the capacitance, V the 
potential difference, and <fi the potential distribution between two conductors, the fol- 
lowing relation is true in general: 

U = | CV 2 = | j* (V0) 2 dr (1) 

The right-hand term describes the energy of an arbitrary electric field bounded by 
the limits of the volume integral and the middle term describes the energy stored in the 
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field of an arbitrary charged capacitor, independent of its geometry. It is understood 
that the volume integral is taken throughout the region bounded by the conductors. The 
symbol e represents the permittivity of the environment (air, vacuum, etc.). It is the 
product of the permittivity of free space and the dielectric constant of the medium between 
the conductors. Equation (1) yields the following: 

c=-%f (V0) 2 d r (2) 

J T 

Thus, the capacitance of an arbitrary condenser charged to a potential difference of 
V volts can be expressed in terms of the potential distribution which exists in the region 
bounded by the plates of the condenser. It is necessary only to find the potential distri- 
bution in the region bounded by the condenser when it is charged. But equation (2) also 
possesses a very remarkable property. Note that it expresses capacitance C as a func- 
tion of potential distribution 0. Any <fi can be inserted in equation (2) to generate an 
expression for capacitance, but obviously only the true 0 will yield the right capaci- 
tance. Of all the possible potential distributions which might exist within a given charged 
condenser, the question is, which potential distribution 0 is actually present? 

From arguments based on the law of the conservation of energy it can be shown that 
if two arbitrary fixed conductors are held at a given potential difference, charges will 
tend to distribute themselves on the surface of the conductors in such a way as to store 
the least amount of energy in the field between them. This is a special case of Thomson's 
theorem (ref. 8). From equation (1) it can be seen that the energy stored in a charged 
condenser is directly proportional to the capacitance of that condenser. It is therefore 
clear that when a fixed potential difference is impressed between two conductors (a con- 
denser) the charges tend to distribute themselves on the conductors in such a way as to 
create a potential distribution corresponding to some minimum capacitance. That is, of 
all the physically realizable values of 0 that might be used in equation (2), the correct 
0 is the one that generates the smallest value for C. This is quite significant.' 

According to the ordinary calculus of extremums, if a continuous function has a minimum, 
then first-order departures from that minimum in the independent variable correspond to 
departures from the minimum in the dependent variable which are only second order to 
the deviations of the independent variable. For example, consider some continuous func- 
tion C(0) which has a minimum at 0 C . Then C(0 c ) is the minimum of that function. 
Now consider a first-order departure from 0 C in 0; that is, consider C(0 c + 60), 
where 60 is some first-order departure from 0 C . From laws developed in the cal- 
culus of extremums it can be shown that C(0 c + 60) differs from C(0 c ) only to sec- 
ond order, even when 60 is a first-order deviation from <p c (ref. 9). Therefore, if 
a first-order approximation to 0 is inserted in equation (2) instead of the true 
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expression for <p, the resulting expression for the capacitance will be a second-order 
approximation to the correct capacitance because, as was determined above, the true 
potential distribution <p corresponds to a minimum capacitance C. 


The preceding analysis has been based in part on arguments developed in the ordi- 
nary calculus of extremums. But <p equation (2) does not represent a variable in the 
ordinary sense; it represents potential distribution, which is itself a function of coordi- 
nate space. It might therefore be argued that an analysis based on the ordinary calculus 
of extremums does not necessarily apply to equation (2). But the analysis can be rigor- 
ously extended to cover the case at hand by invoking the calculus of variations instead of 
the calculus of extremums. This analysis is a bit more involved than the one that has 
been outlined and will not be presented here, but the results are the same as if <fi were 
treated as an ordinary variable instead of a function (ref. 10). 

The important point is that the capacitance of a condenser with a completely arbi- 
trary geometry can be determined from a knowledge of the potential distribution which 
exists when the condenser is charged to some fixed potential difference. Furthermore, 
if the true potential distribution is unknown, a second-order approximation to the capaci- 
tance can be derived by using only a first-order approximation to the true potential 
distribution. 

As a simple example of an application of this analysis, consider the case of a 
parallel-plate capacitor with one plate grounded and the other at a potential V with 
respect to ground. The grounded plate is assumed to be in the z = 0 plane and the 
charged plate is in the z = b plane. The potential distribution, known to be linear, 
might look like the following: 



( 3 ) 


If the plates are rectangular, with length L and width W, equation (2) becomes 
r z=b py=W r*x=L v 2 


r z=a py=w px=j_, v z 
C = — \ \ \ — dx dy dz 

J z=0 J y=0 J x=0 b z 


or 


(5 

This is recognized as the correct formula for the capacitance of a parallel- plate capaci- 
tor with plate dimensions L by W separated by a distance b. 

Now the potential will be assumed to be unknown and the capacitance of a parallel- 
plate condenser will be estimated by inserting a trial potential into equation (2). Of 
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course, the potential distribution of a parallel-plate capacitor is known to be linear, but 
to illustrate the usefulness of equation (2) in deriving close approximations to capacitance 
from rough approximations of potential, a nonlinear potential is inserted into equation (2). 
Assume, for example, that 0 is a sinusoidal potential: 


0 


V sin 


7TZ 

2b 


( 6 ) 


This potential meets the prescribed boundary conditions but is certainly not linear. To 
see what sort of capacitance it produces, equation (6) is inserted into equation (2) to yield 


pz=b r>y=W px=L v 2 2 9 _ z 

C = — \ 1 \ cos 2 — dx dy dz 

y z J z=0 J y=0 J x=0 4b^ 


or 


7 t 2 eLW 
8 b 


Note that the sinusoidal trial potential leads to the correct capacitance formula 
multiplied by a factor slightly greater than 1, namely ir^/s. This is certainly a very 
close approximation to the true capacitance, considering how poor the potential approxi- 
mation was. If other incorrect trial functions are used in equation (2), other approximate 
capacitance formulas will be generated, but just as in this example, the approximation to 
the correct capacitance will always be much closer to the true capacitance than the trial 
potential is to the real potential. Actually, a very accurate expression for capacitance 
can be obtained from a close approximation of the potential; the better the potential 
approximation, the better the expression for capacitance. This method of deriving 
second-order approximations to capacitance from first-order approximations to potential 
is an application of what is called the electrostatic theory of least action. It is a technique 
whereby approximate solutions can be found to otherwise intractable problems. 

In the following sections of this paper, the potential distribution of a grid-plate con- 
denser is determined and used in equation (2) to generate an expression for the capaci- 
tance of a grid-plate condenser. The field solution entails a few approximations and 
therefore is not an exact solution, but because of the variational technique used, the 
resultant capacitance formula is quite accurate. 


Two-Dimensional Grid-Plate Potential Distribution 

The field problem has been cast in a Cartesian coordinate system with the flat plate 
(dimensions L by W) in the z = 0 plane and the grid of wires in the z = b plane. 

The wires have a radius r and are separated a distance a from axis to axis. The 
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plate is assumed to be grounded (i.e., at 0 potential) and the wires of the grid are at 
potential V with respect to the plate. In order to simplify the field problem, the grid 
is imagined to consist of two "subgrids," one subgrid being composed of wires running 
parallel to the X-axis and the other being composed of wires running parallel to the 
Y-axis. The problem is thus simplified into a pair of two-dimensional field problems. 
The solutions to these two problems are combined to yield the potential distribution of 
the three-dimensional grid-plate configuration. Figure 2 illustrates one of the subgrids 



The potential distribution corresponding to figure 2 must satisfy Laplace's equation 
and meet the boundary conditions. The uniqueness theorem guarantees that any such 
solution will be correct and unique. 

Laplace's equation in two dimensions, expressed in Cartesian coordinates, is 

V 2 0(x,y) = 0 (7) 

One possible solution to this equation is 

4> =T> cos Q!X sinh az (8) 

The constants D and a must be evaluated to make equation (8) meet the boundary 
conditions. 

Since the wires are indistinguishable, <p is the same at x = 0 as at x = a for 
any given z. Then 
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cos 0 = cos aa 


aa = 2mr 



<p = D cos 


2n77X 

a 


sinh 


2nnz 

a 


(n - 0, 1,2,3,. . .) 


( 9 ) 

(10) 


Parallel to the X-axis in figure 2, the field changes periodically. This suggests 
that the field can be expressed as a Fourier series. (Indeed, any valid solution to 
Laplace's equation can be expressed as a Fourier series.) Equation (10) has the form 
of the nth term of such a series, which might look like the following: 

OO 

<P = V D n cos 2= sinh *= ( 11 ) 

n=l 


At large distances from the grid, the potential must approach that of a uniform 
charge distribution. That is, in the region between the grid and the plate, as z/a 
approaches zero the potential must approach a linear function of z. This function <fi 0 (z) 
is the potential corresponding to the constant electric field that would exist at large dis- 
tances from the grid. Equation (11) is modified to include the <fi 0 (z) term: 


<p = tp Q { z) + D n cos ^ n ^ rx . sinh 
n=l 


2nTrz 

a 


(12) 


The function <p Q (z) is linear in z , does not depend on x, and is zero when z = 0. 
When z/a approaches zero, the oscillating terms in the sum approach zero and the 
potential approaches 0 O as required. 

Equation (12) still satisfies Laplace's equation. Now D n and 0 O ( z) must be 
determined so that equation (12) will meet the boundary conditions. Before D n and cp 0 
can be found, the boundary conditions must be established. 

The potential is zero at z = 0 by the way the problem is set up. Since 0 o (z=O) = 0 
by definition, equation (12) already satisfies the lower boundary condition. The potential 
on the upper boundary is V wherever there is a wire and something less than V in the 
space between the wires. Figure 3 is a schematic representation of the potential in the 
z = b plane. 
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Figure J. - Potential at upper "boundary. 


Since 0(x, b) is a periodic function, it can be evaluated over the range 0=x=a 
and expanded in a Fourier series. The potential over the range 0 = x = a is the poten- 
tial corresponding to a parallel array of wires of radius r separated center-to-center 
by a distance a with charge residing on one side of the wires (because of the influence 
of the grounded plate). The potential is considered only on the line connecting the centers 
of the wires (the z = b plane). 

The method of images will be used to solve the field of a grounded plate near an 
array of parallel charged wires. The potential of an array of parallel charged wires in 
the presence of a grounded plate is the same as that of an array of parallel line charges, 
each one corresponding to the axial center of charge of one of the wires. If the grounded 
plate is replaced by a parallel array of image line charges, as in figure 4, the field is 
unchanged. For each wire, the location of the axial center of charge will depend on the 
wire radius and on the wire-plate separation. If the distance between each wire and the 
grounded plate is very large, the center of charge can be said to correspond to the axis of 
the wire. If the wire-plate separation is vanishingly small, the axial center of charge for 
each wire will approach the bottom edge of the wire. 

In figure 4, 6 represents the distance that the axis of charge of each wire is dis- 

placed from the geometric axis of the wire due to the influence of the grounded plate. 
Consider the potential at point P(x,b), which lies in the plane of the grid. The potential 
at P(x,b) is composed of contributions from all the line charges and image line charges. 
If d n represents the distance from point P to the nth line charge and d n ' represents 
the distance from point P to the nth image line charge, they can be expressed by 

d n = ((x - na) 2 + 6 2 d n ' = |(x - na) 2 + (2b - 6) 2 
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Figure 4. - Image line charge distribution of a charged grid and grounded plate. 



Both of these expressions hold for positive and negative n. The potential of a line charge 
can be shown to be 

d> = ~-L- In d + Const. 

27 re 


where X is the linear charge density, d is the radial distance from the line charge, 
and the arbitrary constant defines the zero of potential. Therefore the potential at 
point P due to the nth line charge is 

d> n = In d n + Const. 

2ire n 

and the potential at point P due to the nth image line charge is 

0 n ' = In d n ' + Const. 

2776 

The potential at point P is then found by superposing the <p n and $ n ' values: 


cp(x, b) = ■£- ^ jln{(x - na) 2 + (2b - 6) 2 J ^ - lnj(x - na) 2 + S 2 J ^ "1 + Const. 


n=-°° 

+°° 


y 

8tT€ Li 


In ( * - na)2 1 (2b - a > 2 + Const. 


Lu 2 o 

n=-°° (x - na) + 6^ 


(13) 


or, finally, 


^ (x ’ b)= ^i ln TT 


r +°° 

r nl 

j TT 

j , 4b(b - 5) 

1 I 

n=-°° 

(x - na)^ + 5 2 j 


+ Const, 


(14) 


The maximum value of 6 is r, the wire radius, corresponding to the case when 
the plate is almost touching the grid of wires. Of course, the minimum value of 6 is 
zero, when the plate is far away from the grid. If a 2 » r^ (which is true for a large- 
mesh grid of fine wires), 6 can be neglected in equation (14) for all cases except the- 
n = 0 case. Since x can be as small as r (<fi is defined by eq. (14) in the region 
between adjacent wires) and since 6 can be as large as r (corresponding to a vanish- 
ingly small grid-plate separation), it is not immediately obvious that 6 can be neglected 
in equation (14) if n = 0. In appendix A it is shown that even when the plate is close to 
the axis of the wires, neglecting 6 in equation (14) introduces a negligible error in the 
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n = 0 term. The error introduced in the other terms of equation (14) will be even less. 
Therefore, 6 is dropped from equation (14). With 6 assumed to be zero for all 
practical purposes, equation (14) can be restated as follows: 



+ Const. 


This infinite product can be expressed in closed form (ref. 11) as follows: 


(15) 



cosh 2ab - cos ax 
1 - cos ax 


(16) 


where a = — . 

a 


Equation (14) thus becomes: 


0(x,b) = -Uln - c -°g. h 2 .9* - cos + Const. 
8776 \ 1 - cos ax / 


(17) 


To evaluate the constant, consider the fact that when x = r, 0(x,b) must be V, 
the potential of the wire: 

0(r,b) = V = A_(ln cosh 2ab - cos ar) + Const> 

8 ire\ 1 - cos ar / 

Thus 


Const. = V - in cosh 2ab - cos ar 
8776 1 - cos cur 


Then 


0(x,b) 


_X_4 cosh 20b - cos ax j + y _ ln cosh 2ab - cos ar 
8776 V 1 - cos ax / 8776 1 - cos ar 


(18) 


or 


4>(x,b) = V + — ln 
8776 


(cosh 2ab - cos ax)(l - cos ar) 
(cosh 2ab - cos ar)(l - cos ax) 


(19) 
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Of course, equation (19) holds only between the wires since the potential inside a con- 
ductor must be the same as the potential on the surface. The potential inside the wires 
(0 < x < r and (a-r) < x < a) is V. Therefore, for 0 < x < a, $(x,b) is as follows: 






V 

(0=x=r) 

0(x,b) = < 

ir . X (cosh 2Q;b - cos o>x)(l - cos ar) 

(r S x g (a-rj 

87re (cosh 2ab - cos ar)(l - cos ax) 


V 

[(a-r) S x ^ aj 




In equation (12) the potential distribution was cast in the form of a Fourier cosine 
series with two as yet unidentified constants. In order to evaluate these constants, equa- 
tion (20) is expanded in a Fourier cosine series and compared term by term with equa- 
tion (12). The expansion is carried out in appendix B with the following result: 


where 


0(x,b) = V^l -^- + T 


cos ar „ . . ar 

+ 2 ln[2 sin — 


cosh 2o!b 
2 VT[tt - 2 Si(ar)] yT 


cosh 2o?b| 


'cos ax 


+ ; 3YZ n ~ 2 Si(nar) cog nax 
t r n 

n=2 


87reV 


a = 


2-rr 


'X=nar 


px=n 

Si(nar) = \ 

J x=0 


sin x 
x 


dx 


( 21 ) 


Equation (21) represents the grid-plate potential distribution in the plane of the grid. 
It is the upper boundary condition. That is, when equation (12) is evaluated at z = b (the 
plane of the grid) the constants (p Q and D n must be such that equation (12) and equa- 
tion (21) are identical. Note that the upper boundary has been expressed as the sum of a 
<p Q term, a term, and an infinite series of terms <p n , where n = 2. This is simply 
because equation (20) was taken as the sum of three functions, and the Fourier cosine 
expansion of one of them contained only two terms (see appendix C). It is therefore nec- 
essary to express equation (12), the general solution, in a form similar to equation (21) 
in order to evaluate the unknown coefficients. Clearly, equation (12) can be expressed as 
follows: 
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( 22 ) 


0(x,z) = <p o (z) +DiCOS^ sinh ^ + Y D n cos sinh 

n=2 

Consider this expression at z = b: 

OO 

0(x,b) = 0 o (b) + D X cos ^ sinh ^ ^ D n cos sinh 


n=2 


(23) 


This is just the general solution of the potential problem evaluated at the upper boundary. 
If equation (23) is compared term by term with equation (21), the upper boundary condition, 
the following equalities are identified: 


0 o (b) = V<^ 1 - — + T 


cos ar r, , /o „• 

+ 2 In 2 sin — - 

cosh 2ab V 2 


2 7 i x . u 277b jZVT^ir - 2 Si(o-r)] VT 


Di cos sinh o 

i a a 


cosh 2ab | 


>cos ax 


Dn 


cos 


2n7ix 

“a 


sinh 


2n7rb 

a 


2VT tt- 2 Si(nar) 
77 n 


cos nax 


Since 


a 



these equations can be written as 


0 o (z) = v<;i-^ + t 


cos “ r + 2 1n(2sin^ 


cosh 2ab 


z 

b 


(24) 


Di=- 


VT r“ [5 ~ 2 Si(arJ| X 


sinh ab 


cosh 2ab/ 


(25) 


D n = 


2VT 77-2 Si(nar) 


sinh nab n77 


(n ^ 2) (26) 


Equations (24) to (26) are inserted into equation (22) to yield the general solution of 
the potential problem in two dimensions: 
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0(x,z) = V 



+ T 


cos cl r 
cosh 2ab 


+ 2 ln(2 sin ^ 


z 

b 


t 


VT fafft - 2 Si(ar)] 


sinh ceb 


2VT V 77 - 2 Si 

77 L, n 


cosh 2ab 
Si(no!r) sinh naz 


>sinh az cos ax 


n=2 


sinh nab 


cos nax 


( 27 ) 


This is the solution to the problem in two dimensions, but the actual grid consists 
of a plate and two perpendicular arrays of parallel wires. The solution to the three- 
dimensional problem is thus a superposition of a pair of two-dimensional solutions. That 
is, <p(x,y,z) is in reality <p(x, z) (eq. (27)) superimposed with the following equation for 
<p( y,z) (eq. (27) with x replaced by y): 


<My,z) 


v(i + t 


cos ar 


cosh 2ab 
VT fejjr - 2 Si(arj] 


+ 2 ln|2 sin 


ar 


sinh ab 


77 


cosh 2ab| 


>sinh az cos ay 


V 77-2 Si(nar) sinh naz 

+ 2VT > 1 cos nay (28) 

L-j n77 sinh nab 

n=2 


< 


The cp Q term is unaffected by the superposition. It is just the potential which 
exists at large distances from the grid where the oscillating terms in the series do not 
add appreciably to the potential; that is, it corresponds to the uniform charge distribution 
which the grid appears to have at large distances from it. For a given z, <p Q is a func- 
tion only of the charge density and since superimposing two similar (perpendicular) wire 
arrays does not change the overall charge density of the grid, the <p Q term remains 
unaffected. When </>(x,z) and <p( y,z) are superimposed, the result is as follows: 


<Mx,y,z) = V < 1 - + T 


cos ar o i /o • ar 
+ 2 In 2 sin — 


cosh 2ab 


VT ~ ^ Si(ar)] 


sinh ab 


77 


cosh 2abi 


>sinh az(cos ax + cos ay) 


2VT V 77-2 Si(nar) sinh naz, x 

+ / — 1 (cos nax + cos nay) 

v L-i n sinh nab 

n=2 


(29) 
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Equation (29) describes the potential distribution that exists between a flat grounded 
plate and a square-mesh grid of wires charged to some potential V with respect to the 
plate. All the parameters are easily measured except T, which is defined by 


It would be most convenient to be able to express A in terms of some easily mea- 
sured parameters. This can be done as follows. The capacitance of a system consisting 
of a long charged wire parallel to an infinite grounded plate is 


where q is the excess charge on the wire and V the potential difference between the 
wire and the plate. From this it follows that 



where is the capacitance per unit length of the system. Thus the linear charge 
density of a charged wire in the presence of a grounded plate is 


A = VC Z 

The capacitance per unit length of a wire-plate system is (ref. 12) 


Thus 



(30) 


and 


2rreV 



T = 


1 


4 In 


1 

r 



(31) 


(32) 
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Equation (32) can be inserted into equation (29) to find the potential distribution of a 
charged grid-plate capacitor in terms of easily measured parameters, namely the grid- 
plate spacing b, the parallel -wire separation a, the wire radius r, and the grid-plate 
potential difference V. 


Computation of Capacitance From Potential 

The fact that the capacitance between arbitrarily shaped conductors can be derived 
from the potential distribution that exists when the conductors are charged was discussed 
in the section entitled "Approach." The relation between capacitance and potential was 
expressed in equation (2): 



This volume integral is taken throughout the region bounded by the conductors. An 
expression for <fi is derived in the previous section (eq. (29)), and in appendix C this 
expression is inserted into equation (2) to yield 


where 




2 nab sinh 2nab\ 


n=l 




r 


cos ar 0 . /„ . ar 

+ 2 In 2 sm 


cosh 2ab 


(33) 


Al = 


T fefn - 2 Si(ar)] 


sinh ab 


cosh 2ab I 


t 


An = 


2T 77-2 Si(nar) 
sinh nab htt 



T = 


1 

4 In 

/ \ 

•iNb^ - r^ + bj 


( n = 2) 
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Note that in equation (33) the term outside the parentheses is the formula for the 
capacitance of a parallel solid-plate capacitor with plate dimensions L and W and 
separation distance b. The term in parentheses is a measure of how this capacitance 
is changed when one of the plates is replaced by a square-mesh grid of wires with radius 
r separated axis-to-axis a distance a. 

It should be remembered that the expression used for the potential distribution to 
determine the capacitance was only an approximate one. Although the capacitance formula 
is derived by a variational method that requires only an approximate expression for <p to 
yield a very close approximation to the true capacitance, nevertheless equation (33) repre- 
sents only an approximation to the true grid-plate capacitance. Just how good this approx- 
imation is can be determined by actually measuring the capacitance of a number of grid- 
plate configurations and comparing the results with predictions based on equation (33). 
Values reported in the next section reflect equation (33) carried out (by computing calcu- 
lator) to 15 terms. Actually, the series in equation (33) converges quite rapidly, and very 
good results can be had by retaining fewer terms in the series. 

EXPERIMENTAL CONFIRMATION OF THEORY 

Capacitance measurements were made on a number of grid-plate configurations. 

The grid-plate capacitors were composed of a single flat polished plate of aluminum and 
a number of interchangeable grids. Each grid in turn was positioned above the aluminum 
plate by means of a number of columns of small cardboard spacers, which varied in height 
from 5.1 mm to 12.7 mm in 1.3-mm increments. Several data points were taken for each 
grid, every point corresponding to a different grid-plate separation. Thus, for each grid, 
an experimental curve of capacitance as a function of grid-plate separation was generated. 
A theoretical curve was obtained by programing the capacitance formula on a computer 
with instructions to increment b, the grid-plate separation parameter. Each computer- 
generated data point was also plotted. Comparisons between the empirical and theoretical 
curves indicate that the grid-plate capacitance formula is exact within experimental error. 

Figure 5 illustrates a typical experimental configuration. The flat plate is mounted 
on a hollow block of polystyrene to displace it from the table and reduce the effect of stray 
capacitance. The spacers are positioned on the plate and the grid is placed on top of the 
spacers. The spacers were positioned so as to prevent the grid from sagging. End 
effects were suppressed by using relatively small grid-plate separation distances and 
relatively large grids and plates (ratios of plate area to separation on the order of 500 
to 1). The capacitance measurements were made with a precision industrial capacitance 
bridge. The recorded data, corrected for effects of the cardboard spacers, fringing 
fields, and meter -lead capacitance, are reported in tables I and H for two grid-plate 
capacitors. In each case, predicted and measured values are given for a number of grid- 
plate separations. 
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Figure 5. - Configuration for capacitance measurements. 


TABLE I.- 6.4-mm MESH 


TABLE H.- 25.4-mm MESH 


t 


Wire radius, 0.318 mm ± 0.013 
Mesh dimension, 6.350 mm ± 0.025 
Edge of grid and plate, 30.5 cm ±0.1 
Spacer thickness (each), 1.270 mm ± 0.010 


Wire radius, 0.635 mm ± 0.051 
Mesh dimension, 25.400 mm ± 0.102 
Edge of grid and plate, 30.5 cm ±0.1 
Spacer thickness (each), 1.270 mm ± 0.010 



Grid-plate capacitance, 

Number of 

picofarads 

spacers 

Experimental 

Theory 

4 

139.8 ± 1.5 

142.8 ± 3.5 

5 

116.3 ± 1.9 

116.4 ± 0.3 

6 

99.7 ± 0.8 

100.3 ± 0.7 

7 

88.2 ± 1.0 

88.7 ± 1.0 

8 

79.0 ± 0.6 

79.1 ± 0.3 

10 

66.0 ± 0.6 

65.3 ± 0.5 



Grid-plate capacitance, 

Number of 

picofarads 

spacers 

Experimental 

Theory 

4 

108.3 ± 6.5 

102.8 ± 1.3 

5 

95.3 ± 2.0 

90.2 ± 1.0 

6 

83.9 ± 2.1 

81.1 ± 0.9 

7 

74.7 * 1.5 

74.0 ± 0.7 

8 

68.0 ± 1.1 

68.4 ± 0.8 

9 

61.4 ± 0.8 

63.8 ± 0.7 

10 

58.0 ± 0.8 

60.0 ± 0.7 
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Figures 6 and 7 are computer-generated curves of the capacitance formula, equa- 
tion (33), as applied to the grid-plate capacitors described in tables I and II. Superim- 
posed on the theoretical curves are measured data points, with error bars. 

DISCUSSION OF RESULTS 

All the uncertainties in the experimental values reported, both in geometric mea- 
surements and measurements of capacitance, are the standard deviations based on 10 mea- 
surements per point. Uncertainties in the theory reflect the standard deviations in mea- 
sured grid-plate geometry parameters. 

In the theoretical analysis, the thin-wire approximation was made a number of 
times. Hence, the analysis can be accurately applied only to geometries in which the 
radius of the grid wires is small compared with other dimensions such as the grid-plate 
spacing and the distance between adjacent parallel wires. 

The charge distribution on each wire is slightly affected by the presence of charge 
on neighboring wires, if the grid is assumed to be finite. However, the effect is quite 
small everywhere except near the very edges of the condenser, and even there the effect 
is small for thin wires, since the center of charge of each wire is constrained to lie some- 
where inside the wire. 

The axial center of charge of each wire in the grid was assumed to coincide with the 
geometric axis of each wire. Because of the influence of the grounded plate, this is not 
really the case; the charge is shifted slightly toward the plate. However, it is shown in 
appendix A that even when the plate is very close to the grid, the shift in the axial center 
of charge is quite small and can in fact be neglected for a wide range of grid-plate 
separations. 

In making the capacitance measurements, the grid of each grid-plate condenser was 
suspended over the plate by means of cardboard spacers. Care was taken in the place- 
ment of the spacers to minimize the effects of sag in the grid; however, the data may be 
biased slightly high where the grid-plate spacing is relatively small and grid sag would 
have a larger percentage effect. Care was also taken to minimize the effects of stray 
capacitance by isolating the grid-plate condensers as much as possible; however, some 
small stray capacitance effects are inevitable. They would also tend to bias the data 
slightly high. Finally, the grids used in the measurements were made out of metal which 
had been dipped in a galvanizing solution. Consequently, there were some irregularities 
in the grid, especially at the wire intersections where the galvanizing solution tended to 
concentrate. In measuring the wire radii, several measurements were made with a 
micrometer at different places on the wire and the numbers listed for wire radii are 
average values. 
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These errors are all small, and care has been taken to minimize them. In any 
case, some of the errors tend tc cancel so that the cumulative effect of all the error 
sources is very small. 

The theory predicts a slightly modified inverse-first-power dependence of grid- 
plate capacitance on grid-plate separation, which is verified by the experimental data. 

This result is not surprising, since the parallel grid-plate capacitor is in effect a paral- 
lel solid-plate capacitor with one plate possessing a diminished surface area. Even 
though the shape of the curve of capacitance as a function of separation is not surprising, 
the magnitudes of the capacitance values predicted (and measured) are quite unexpected. 
This is because the relative surface area of a fine-wire open-mesh grid is essentially 
negligible compared with the surface area of a solid plate that has the same length and 
width. Since the capacitance of a parallel-plate capacitor is directly proportional to the 
surface area of the plates, replacing one of the plates with an open-mesh grid of wires 
having a relatively minute surface area would be expected to cause a drastic reduction in 
the capacitance. Such was not found to be the case. Table in compares measured capac- 
itance values from tables I and n with a parallel solid-plate condenser of the same length 
and width. Figure 8 is a graph of these data. 


TABLE HI.- CONDENSER COMPARISONS 


Number of 

Capacitance, picofarads 

1.3-mm 

25.4-mm 

6.4-mm 

Parallel 

spacers 

mesh 

mesh 

plate 

4 

108.3 

139.8 

188.6 

5 

95.3 

116.3 

143.8 

6 

83.9 

99.7 

120.4 

7 

74.7 

88.2 

103.2 

8 

68.0 

79.0 

92.0 

9 

61.4 

71.2 

82.0 

10 

58.0 

66.0 

74.5 


These data illustrate that the capacitance between a plate and an open-mesh grid is 
surprisingly similar in magnitude to the capacitance of a parallel solid-plate condenser 
even when the surface area of the mesh is less by orders of magnitude than the plate area. 
This has led to the conclusion that the minute surface areas which characterize such open- 
mesh structures as the LOFT radio telescope (described in the introduction) do not impose 
a serious restriction on the use of capacitive, noncontacting displacement-measuring 
transducers in the modal analysis of the dynamic response of such structures. 
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POTENTIAL APPLICATION 


A laboratory prototype capacitive transducer system has been developed for appli- 
cation to the LOFT structural-stability tests mentioned in the introduction. In laboratory 
prototype testing, a 15-cm-diameter circular probe mounted 30 cm from the LOFT reflec- 
tor net has detected displacements of the net in the range from 0.05 cm to 20 cm. Deflec- 
tions in the LOFT reflector net on the order of 0.2 cm have been detected from as far 
away as 75 cm. Development of calibration techniques and further refinements would be 
required prior to application of this instrument system. Other potential applications 
include vibratory and modal analysis of small-surfaced hot structures and displacement 
measurements on fragile or lightweight surfaces in general. 

CONCLUDING REMARKS 

The purpose of this paper has been to demonstrate the feasibility of employing a 
capacitive displacement-measuring transducer in the modal analysis of open-mesh grid 
structures possessing minute surface areas. The analysis has been centered about a 
mathematical model of a square-mesh parallel grid-plate capacitor. An expression for 
the capacitance of such a capacitor was derived, and values based on this expression 
were verified by experiment. It has been concluded that the capacitance between a flat 
plate and an open-mesh grid of fine wires is of sufficient magnitude to justify the develop- 
ment of a capacitive transducer for application to displacement-measurement problems 
involving open-mesh structures. A program has been initiated to develop such a trans- 
ducer, and a laboratory prototype system has been built. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 23, 1971. 
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APPENDIX A 


DETERMINATION OF THE ORDER OF MAGNITUDE OF 6 

By Norman C. Wenger 
Lewis Research Center 

The axis of charge of a long, uniformly charged wire coincides with the axis of the 
wire. If an infinite grounded plate is brought near a charged wire, the axis of charge will 
be shifted somewhat by a radial distance 6, where 6 is greater than zero and less than 
the radius of the wire. 

The center of the wire is assumed to be located a distance z = b above a grounded 
plane and the axis of charge is shifted a distance 6 below this. Then 6 is calculated 
by replacing the grounded plane with an image line charge (see sketch (a)): 

0(P) = In di + — - In do + Const. 

2-ne 2-ne 


Z 



Since dj = d2 at z = 0, it follows that <£( z=0) = 0 for all x, and hence the constant 
equals zero. Then 
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APPENDIX A - Continued 


tf>(P) = 


X 

27 re 



When d 2 /di is a constant, 0(P) is a constant. The (equipotential) surface that corre- 
do 

sponds to — = k (a constant) is determined from the following derivation: 

<*1 


dj 2 = d 2 + 6 2 - 2d6 cos^0 

d 2 2 = d 2 + (2b - 6) 2 - 2d(2b - 6)cos^0 + 
d x 2 = d 2 + 6 2 - 2d6 sin 9 

d 2 2 = d 2 + (2b - 6) 2 - 2d(2b - 6)sin 6 

.2 m 2 

d 2 = kdj 

d 2 + (2b - 6) 2 + 2d(2b - 6)sin 6 = k 2 d 2 + k 2 6 2 + 2d6k 2 sin 6 
(1 - k 2 )d 2 + (2b - 6) 2 - k 2 6 2 + (4b - 26 - 26k 2 )d sin 9 = 0 


Any value of k corresponds to an equipotential surface. If 4b - 26 - 26k 2 = 0 
or k 2 = 2 ^ " this surface is a cylinder of radius 


r = 


I (2b - 6) 2 - k 2 6 2 
k 2 -l 


centered at z = b. After substitution of k 2 = 

6 


r 2 _ (2b - 6) z - (2b - 6)6 _ 


= 2b 6 - 6^ 


2b - 6 


- 1 


Therefore 


6 2 - 2b6 + r 2 = 0 
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APPENDIX A - Concluded 


Solving for 6 gives 


6 = b ± \J 

The minus sign is chosen: 



The smaller the ratio b/r, the larger 6 will be. The smallest value of b/r 
measured in this paper corresponds to 6 = 0.06r. In the body of the paper it was nec- 
essary to compare r^ and 8% to see whether 5^ was negligible when compared with 


62 


r In the worse case, — < 0.004. Clearly, then, 


r" r 

(See eq. (14) and the discussion immediately following it.) 


1 + V‘ 


is a good approximation. 
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APPENDIX B 


FOURIER COSINE SERIES EXPANSION OF THE TWO-DIMENSIONAL 
UPPER BOUNDARY CONDITION 


In this appendix, $(x,b) as derived in equation (20) will be expanded in a Fourier 
cosine series. Equation (20) can be written 


fv 


<Mx, b) 


J v + VT ln (cosh 2a b - cos ax)(l - cos ar) 
] (cosh 2ab - cos ar)(l - cos ax) 


l v 


The function $(x, b) has the following form: 


V 

(p{x, b) = ^ V + VTF(x) 

V 

V. 

where 


(0 ^ x S r) 

[r = x = (a-rj (Bl) 
[(a-r) = x = a] 

(0 = x = r) 

[r = x = (a-r)] (B2) 

[(a-r) = x = a] 


F(x) = ln 


(cosh 2ab - cos ax)(l - cos ar) 
(cosh 2ab - cos ar)(l - cos ax) 


(B3) 


A straightforward evaluation of the Fourier coefficients for the function cp(x,b) 
involves the solution of some integrals that cannot be evaluated in closed form. In order 
to avoid this difficulty, <p(x, b) has been synthesized from three functions which are 
easier to treat by Fourier analysis. The first function is called M(x) and is defined as 
follows: 


M(x) = V [l + TF(x] 

where F(x) is defined by equation (B3). 
|M(x) 


(0 = x = a) 

The second function is G(x): 

(0 § x g r) 


G(x) = < 0 

M(x) 


jr = x = (a-r)] 
[(a-r) = x S a] 


(B4) 


(B5) 
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APPENDIX B — Continued 


The third function is H(x), a simple step function: 


H(x) = 


(v 



(0 = x = r) 

jr = x = (a-r)j (B6) 
j(a-r) = x = a] 


Figure B1 illustrates how M(x), G(x), and H(x) combine to form ©(x,b): 


<p(x,b) = M(x) - G(x) + H(x) 


(B7) 


The Fourier coefficients of M(x), G(x), and H(x) will be derived and combined to form 
a Fourier series representation of cf>(x,b). 

Begin by considering the function F(x), defined by equation (B3), which can be 
expressed as follows: 


F(x) = In 


cosh 2ab - cos ax 
cosh 2o;b - cos ar 


+ In 


1 - cos ar 
1 - cos ax 


(B8) 


The first term on the right of equation (B8) can be expressed as follows: 

cos ax 


, cosh 2ab - cos ax , 

In = In 

cosh 2ab - cos ar 


1 - 


cosh 2ab 
cos ar 
cosh 2ab 


In 


1 - 


cos ax 
cosh 2ab 


1 + 


cos ar 
cosh 2ab 


In 1 + 


cos ar - cos ax 
cosh 2ab 


(B9) 


or 


cosh 2ab - cos ax _ cos ar cos ax 
cosh 2ab - cos ar cosh 2ab cosh 2ab 


(BIO) 


The first term on the right of equation (B8) has been reduced to a two-term Fourier 
cosine series by making successive first-order expansions of terms of the form 1 ± v, 
where v « 1. 
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APPENDIX B - Continued 


The second term on the right of equation (B8) can be expressed as follows: 

1 „ 1 - cos ar = 1 In f 1 - COB ar f = i ln(l - cos ar)2 - - ln(cos ax - l) 2 (Bll) 

1 - cos ax 2 \1 - cos axj 2 2 

The following equality comes from page 358 of reference 13: 



(B12) 


and can be expressed as 

OO 

4 V cos nax _ ln 4 

Lj n 

n=l 

OO 

2 £ £OS_nax . ln 2 ( B13) 

n=l n 

Equation (B13) is inserted into equation (Bll) to obtain 

OO 

t 1 - cos ar o V cos nax , 0 1 , N 2 

ln = 2 / + in 2 + - ln(l - cos ar) (B14) 

1 - cos ax Ay n 2, 

n=l 


ln[(cos ax) - f] = - 


or 


— ln(l - cos ax)^ 
2 


Equations (B14) and (BIO) are inserted into (B8) to yield, after some algebraic 
manipulation, 


F(x) 


cos ar , L . o ar\ 

= + In 4 sin z — 

cosh 2ab y 2 j 


OO 

cos ax g V cos no?x 

cosh 2ab A, n 
n=l 


(B15) 


This function is then used in equation (B4) to form the Fourier series representation of 
M(x): 


M(x) = V + VT 


cos ar 
cosh 2ab 


+ ln (4 sin^ 



OO 

VT cos ax + 2VT V cos nax 

cosh 2ab A. n 

n=l 


(B16) 


The next function to be analyzed is 
where G(x) is nonzero, 


G(x), defined in equation (B5). 


In the region 
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APPENDIX B - Continued 


cosh 2ab - cos ax ~ cosh 2ab - cos ar 


Therefore F(x), defined in equation (B3), can be expressed to a good approximation as 


F(x) ~ In 


1 - cos ar 
1 - cos ax 


[b < x < r, (a-r) < x < a| 


From equations (B17), (B4), and (B5), 


V + VT In 


G(x) * {0 


V + VT In 




1 - cos ar 
1 - cos ax 


1 - cos ar 
1 - cos ax 


The general form of the cosine expansion of G(x) is: 

OO 

G(x) = — + ^ a n cos nax 


n=l 


where 


“o - ! £ G < x ) dx 

2 r a 

a n = — \ G(x)cos nax dx 
a J n 


(0 < x < r) 
(r < x < (a-r)) 
[(a-r) < x < a] 


(B17) 


(B18) 


(B19) 


First a Q will be evaluated. Figure B1 shows that G(x) is a symmetric function. 
It is symmetric about a/2; that is, G^| + xj = G^ - xj for 0 < x < a. This fact can be 

verified by direct substitution into equation (B18). 

Since G(x) = 0 for r < x < (a - r), ao can be evaluated as follows: 


l ° " i I 


r V + VT In 1 - cos ar dx + ° 


1 - cos ax 


f£ 


a-r 


V + VT In 1 - cos ar dx 
1 - cos ax 


(B20) 
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APPENDIX B — Continued 


Because G(x) is symmetric about a/2, equation (B20) can be written as 


a Jq 


V + VT In 


1 - cos ar 
1 - cos ax 


dx 


(B21) 


Since a = — and because of the thin-wire approximation r « a, equation (B21) 

ci 

can be put into a more tractable form. Since 

e 2 

cos 6 ~ 1 

2 

for small 8, 


ln 1 - c ° s ° r * ln /g£f = 2 m I 

1 - cos ax \otx 


x 


Thus, equation (B21) becomes 
4 / 




V + 2VT In t- dx 


or 


ao = ~~(1 + 2T) 

a 


(B22) 


Similar arguments apply to the evaluation of a n (i.e., symmetric G(x) and thin- 
wire approximation). Therefore, 




V + 2VT In — jcos n«x dx 


(B23) 


After evaluating this integral, the result is 


a n = -^^Isin nar + 2T Si(no!rj] 
2n7T ^ -* 


(B24) 


where 


rli=na r . 

Si(nar) s f dp 

J M= 0 M 


is the sine integral function, which is tabulated in the literature (see, for example, ref. 14). 


36 



APPENDIX B — Concluded 


Since equations (B22) and (B24) are the Fourier coefficients of G(x), the Fourier 
series of G(x) is 


OO 

G(x) = -?XE(l + 2T) = V -^(sin(no!r) + 2T Si(no!r)|cos ndx (B25) 

a ^ 2n7r 

n=l 


The final function to be analyzed is H(x), defined in equation (B6). For the cosine 
expansion of H(x), the Fourier coefficients are derived in the same way they were 
derived for G(x). They are 


ao - 


4Vr 

a 


an = 


2V sin no?r 
nir 


Therefore, 

CO 

H(x) = + £ 2V _ |!L Sg£ cos nax (B26) 

n=l 


The Fourier series of M(x), G(x), and H(x) have now been derived in equa- 
tions (B16), (B25), and (B26), respectively. When these equations are inserted into equa- 
tion (B7), the result is a cosine expansion of <£(x,b). Then, after some algebraic 
manipulation, 


<Mx, b) = V 



4Tr 

a 


+ T 


cos ar 
cosh 2o!b 



2VT [w - 2 Si(o!rJ 


cos ox + ^ 2V T * - 2 Si(nar) CQg nffX (B27) 

n=2 

This series was programed into a computing calculator and plots of <p{x, b) as a 
function of x were made for a variety of geometric parameters, keeping 15 terms of the 
series. These plots coincided with plots of <p{x,b) as it is expressed in equation (Bl), 
verifying that equation (B27) is indeed the Fourier cosine expansion of equation (Bl). 
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APPENDIX C 


DERIVATION OF GRID-PLATE CAPACITANCE FROM 
THE GRID-PLATE POTENTIAL DISTRIBUTION 

Equation (2), repeated here for convenience, relates the capacitance of an arbitrary 
condenser to the potential distribution inside that condenser: 


C = 


y2 xJ T 


J (*4>y 


'dT 


For a parallel grid-plate condenser, the potential distribution is described by 
equation (29): 


<p{x,y,z) = V' 


r 


- 

) / 

1 - 4Tr + T 

cos ar 2 J /, 

sin — 

\ z + VT J 

a 

cosh 2ah \ 

2 / 

( b sinh ab ] 


— 

— * 

1 1 


cosh 2ab| 


■sinh o!z(cos ax + cos ay) 


w 

2VT V 77-2 Si(nar) sinh naz/„„„ N 

/ * L (cos nax + cos nay) 

T! t- 1 n sinh nab 


n=2 


It is necessary to insert equation (29) into (2) to find the capacitance of the grid-plate 
condenser. 


Since 0(x,y,z) is an infinite series of the form 

OO 

0 - $o + $1 + ^n 


(Cl) 


then 


n=2 


V0 = V0 O + V0J + ^ V0 n 
n=2 


and 


(V0) 2 = (V0 O ) 2 + (V0i) 2 + ( £ V0 n j + 2Vc)) 0 • V0! + 2V0 O • £ V<p n + 2Vc#> 1 • £ V0 n 


\n=2 


n=2 


n=2 


(C2) 
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APPENDIX C — Continued 


The third term on the right side of equation (C2) can be more clearly expressed as 
the sum of two series: 



2 ( v 4>n ) 2 +J^V0 n V0 m 

n=2 n m 

n#m 


Equation (C2) thus becomes 


(C3) 


oo oo oo 

(V0) 2 = (V<£ 0 ) 2 + (V0!) 2 + £ (V<p n ) 2 + 2 V0 O • V0J + 2 V<p o • ^ V0 n + 2V<^ 1 • ^ V0 n 

n=2 n=2 n=2 

OO OO 

+ V ^ V( ^m < C4 ) 

n m 
n^m 


Since each of the terms on the right side of equation (C4) is composed of products 
of complete functions, if the grid is assumed to be composed of an integral number of 
meshes (a condition which insures the conservation of integral periodicity), the integrals 
of the fourth, fifth, sixth, and seventh terms on the right side of equation (Cl) are zero by 
orthogonality. What remains is 



(V<£) 2 dr = 



2 (V<£ n ) 2 dT 
n=0 


(C 5) 


where 


<Po 


A 0 Vz 

b 


(C6a) 


c p n = A n V sinh no!z(cos nax + cos nay) 


(C6b) 


with 


■^o 


4rT T f cos ar 2 
a cosh 2ab 



Al = 


T ~ 2 Si(ar)] 


sinh ab I 


cosh 2 abl 
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APPENDIX C - Continued 


A n - 


2T 7 T - 2 Sj(nttr) 
sinh nab mr 



(nS2) 


T = 


1 


4 In 


1 

r 



From equations (C6), 



and 


(C7a) 


2 2 2 n 

(V0 n ) = (naVD n sinh naz; (sin^nax + sin^nay) + (no!VD n cosh naz) (cos nax + cos nay) 2 


(C7b) 

The plate and grid are assumed to have a length L and a width W and to be 
spaced a distance b apart. Therefore L, W, and b represent the limits of integra- 
tion in the x-, y-, and z-direction, respectively. Equation (C5) becomes 


f (V„V= Y f*~ f =y f LM (V^dxdydz 

J r L -' r . J z-0 J y=0 J x=0 

n=0 J 

Equations (C7a) and (C7b) must be integrated: 

pb=z r W=y pL=x 9 V 2 A n 2 LW 

(V0 o ) 2 dx dy dz = 2 


; z=0 w y=0 ^x=0 


pu=z nw=y pij=x 9 V 2 LWna!D 2 sinh 2nab 

(V0 n ) 2 dx dy dz = ^ 

^z=0 ^v=0 ^x=0 ^ 


->b=z pW=y pL=x 

'y 

Combining equations (C8) and (C9) yields 

/ 2 00 


y <v«) 


2 dr = V 2 LW 


^ + y 

i b „4 


naA n sinh 2nab 


(C8) 


(C9a) 


(C9b) 


(CIO) 
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APPENDIX C — Concluded 


Equation (C 10) is inserted into equation (2) to obtain 


r _ LWe 
b 



2 nab sinh 2no!b| 


v n 


(Cll) 


The term outside the parentheses is the well-known formula for the capacitance of 
a parallel- plate condenser. The expression inside the parentheses represents the change 
in capacitance when one of the plates is replaced by a square-mesh grid of wires of 
radius r, separated by an axis-to-axis distance a. The series converges very rapidly 
and only a small number of terms need be retained to insure a high degree of accuracy. 
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